Extending cosmology: the metric approach 
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Introducing a fundamental constant of nature with dimensions of acceleration into the theory 
of gravity makes it possible to extend gravity in a very consistent manner. In the non-relativistic 
regime a MOND-like theory with a modification in the force sector is obtained, which is the weak- 
field limit of a very general metric relativistic theory of gravity. Since the mass and length scales 
involved in the dynamics of the whole universe require small accelerations of the order of Milgrom's 
acceleration constantao, it turns out that the relativistic theory of gravity can be used to explain 
the expansion of the universe. In this work it is explained how to built that relativistic theory of 
gravity in such a way that the overall large-scale dynamics of the universe can be treated in a pure 
metric approach without the need to introduce dark matter and/or dark energy components. 
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I. INTRODUCTION 

Current cosmological data is generally explained in- 
troducing two unknown mysterious dark components, 
namely dark matter and dark energy. These ad hoc 
hypothesis represent a big cosmological paradigm, since 
they arise due to the fact that Einstein's field equations 
are forced to remain unchanged under certain observed 
astrophysical anomalies. 

A natural alternative scenario would be to see whether 
viable cosmological solutions can be found if dark un- 
known entities are assumed non-existent. The price to 
pay with this assumption is that the field equations of 
the theory of gravity need to be extended and so, new 
Friedmann-like equations will arise. The most natural 
approach to extend gravity arises when a metric exten- 
sion /(i?) is introduced into the theory [see e.g. [t| and 
references therein]. 

In a series of recent articles, Bernal et al. [1, Q , Car- 
ranza, Mendoza, and Torres [lfj, Hernandez, Jimenez, 
and Allen [ljj], Hernandez et alfisj. Mendoza et al. [22 1 
have shown how relevant the introduction of a new fun- 
damental physical constant arj ~ 10 _10 m/s 2 with dimen- 
sions of acceleration is in excellent agreement with differ- 
ent phenomenology at many different astrophysical mass 
and length dimensions, from solar-system to extragalac- 
tic and cosmological scales. The introduction of the so 
called Milgrom's acceleration constant ao in a description 
of gravity means that any gravitational field produced by 
a certain distribution of mass (and hence energy) needs 
to incorporate the acceleration ao together with Newton's 
gravitational constant G and the speed of light c in the 
description of gravity. 



In section |TT] it is shown, through a description of an 
extended Newtonian gravity scenario, the advantages of 
working with a modification of gravity dependent on the 
mass and lengths associated with the dimensions and 
masses of the sources that generate the gravitational 
field, and not with the dynamical acceleration they pro- 
duce on test particles. Section IIIII describes how it is 
possible to build a metric theory of gravity which gen- 
eralises the extended Newtonian description mentioned 
in section [TT] and section IIVI interconnects this extended 
relativistic description of gravity with a metric descrip- 
tion of gravity for which the energy-momentum tensor 
appears in the gravitational field's action. On section IVl 
we use the developed theory of gravity for cosmological 
applications in a dust universe and see how it is a co- 
herent representation of gravity at cosmological scales. 
Finally on section I VII we discuss the consequences of the 
developed approach of gravity and some of the future 
developments of the theory. 



II. EXTENDED NEWTONIAN GRAVITY 

Milgrom [H, [H, [HI constructed a MOdified Newto- 
nian Dynamics (MOND) theory, based on the introduc- 
tion of a fundamental constant of nature ao = 1.2 x 
10~ 10 ms~ 2 in such a way that the acceleration experi- 
enced by a test particle on a gravitational field produced 
by a point mass source M is such that: 



GM 



VapGM 



for 
for 



a > a , 
a <C a , 



(1) 
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where r is the radial distance to the central mass. In 
other words, for accelerations a> ao, Newtonian gravity 
is recovered and new MONDian effects are expected to 
appear for accelerations a < ao- The strong s « ao 
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MONDian regime means that Kepler's third law is not 
valid since for a circular orbit about the central mass 
M, the acceleration a = v/r, where v is velocity of the 
test mass, and so v — (apGM) 1 ^ 4 oc Af 1 / 4 , which is the 
Tully-Fisher relation [sec e.g. 3.2) for the case of a spiral 
galaxy and is the same relation experienced by wide-open 
binaries [17| and by the tail of the "rotation curve" in 
globular clusters [HI [l6[ ■ 

In order to interpolate from the strong a > ao New- 
tonian regime to the weak a <C ao one, the traditional 
MONDian approach is to construct a somewhat built- 
by-hand interpolation function fi(y) in such a way that 



where 



Ms/) 



afj,(y) 



for 
for 



GM 



(2) 



2/> 1, 



and 



a 

V ■= — • 
ao 



The usual approach to MOND as expressed by equa- 
tion @ means that Newton's 2nd law of mechanics needs 
to be modified [see e.g. [lj. As explained by Mendoza 
et al. [22| , a better physical approach can be constructed 
if the modification is made in the force (gravitational) 
sector. Indeed, by the use of Buckingham's theorem of 
dimensional analysis [cf. HH, the gravitational accelera- 
tion experienced by a test particle is given by 



a = a g(x), 

where the dimensionless quantity 



l_M 

r 



and a mass-length scale 



Im 



GM\ 1/2 
ao / 



(3) 



(4) 



(5) 



The length Im plays an important role in the de- 
scription of the theory and is such that when Im 3> r, 
the strong Newtonian regime of gravity is recovered and 
when Im <C r the weak MONDian regime of gravity ap- 
pears. As such, the dimensionless acceleration (or tran- 
sition function) g{x) is such that: 



whose length is much smaller than its mass-length scale 
is in the Newtonian regime (since x 1). The case 
x = 1 can roughly be thought of as the point where the 
transition from the Newtonian to the MONDian regime 
occurs. 

A general transition function g(x) was built by [22j tak- 
ing Taylor expansion series about the correct MONDian 
and Newtonian limits, yielding: 



l±x n 



l±x n 



(7) 



This non-singular function converges to the correct ex- 
pected limits of equation ([6]) for any value of the param- 
eter n > 0. As shown in Figure [TJ the transition function 
g(x) rapidly converges to the limit "step function" 



0O) 



for < x < 1, 
for x > 1, 



(8) 



when n > 3. The parameter n needs to be found em- 
pirically by astronomical observations. The value found 
by Mendoza et al. [22| for the rotation curve of our 
galaxy is n > 3 and the one found by Hernandez and 
Jimenez [l5j . Hernandez, Jimenez, and Allen (l6| . Her- 
nandez, Jimenez, and Allen [l7j is n > 8, with a mi- 
nus sign selection on the numerator and denominator 
on the right hand side of equation ([7]). These authors 
have shown that a large value of n is coherent with so- 
lar system motion of planets, rotation curves of spiral 
galaxies, equilibrium relations of dwarf spheroidal galax- 
ies and their correspondent relations in globular clusters, 
the Faber- Jackson relation and the fundamental plane of 
elliptical galaxies as well as with the orbits of wide binary 
stars. The n = 3 model in which a small, but measurable 
transition is obtained, has also been tested on earth and 
moon-like experiments by Meyer et al. [24[ and Exirifard 
[l2| respectively, showing that it is coherent with such 
precise measurements. In fact, these experiments also 
validate all n > 3 models. 

Care must be taken when the introduction of a new 
fundamental constant of nature with dimensions of ac- 
celeration ao is made. In fact, the introduction of ao 
does not impose any kind causality arguments such as 
the ones given by the velocity of light c. In fact, one may 
think of ao as a fundamental constant needed to transit 
from one gravity regime to another. In this respect for 
example, instead of using ao as a fundamental constant, 
one may define 



a fx, when x » 1, 

— = g(x) := < (6) 
ao I x, when i<1. 

In general terms, a mass distribution whose length is 
much greater than its associated mass-length Im is in the 
MONDian regime (since x 1) and a mass distribution 



So :=a /G= 1.8kgm- 2 , (9) 

as the new fundamental constant of nature. The constant 
S , with dimensions of surface mass density, enters in the 
description of the gravitational theory in such a way that 
equations ^ and are given by: 
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fundamental constant to use since it is directly measured 
through the flattened rotation curves of spiral galaxies. 

The extended Newtonian model of gravity presented in 
this Section is equivalent with MOND on spherical and 
cylindrical symmetry but deviates considerable from it 
for systems away from this symmetry [22[ • As we have 
already shown, there are however many advantages of 
using this approach, the most objective meaning that 
the modification is made on the force sector and not a 
modification on the dynamics. 
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FIG. 1. The figure taken from Mendoza et al. [22| shows 
the acceleration function a in units of Milgrom's acceleration 
constant ao as a function of the parameter x. The thick dash- 
dot curve is the extreme limiting value n — > oo, i.e. a/ao = 
x for x < 1 and a/ao = x 2 for x > 1. The curves above 
and below this extreme acceleration line represent values of 
n — 4, 3, 2, 1, for the minus and plus signs of equation (0) 
respectively. The extreme limiting curve has a kink at x = 1. 



III. RELATIVISTIC METRIC EXTENSION 

Finding a relativistic theory of gravity for which one 
of its non-relativistic limits converges to MOND yields 
usually stran ge a ssumptions and/or complicated ideas 
[see e.g. |3j, [a [29(. A good first approach was provided 
by a slight modification of Einstein's field equations by 
Sobouti (35|, but the attempt is not complete. 

In order to find an elegant and simple theory of gravity 
for which a MONDian solution is found, Bernal et al. 
Q used a metric correct dimensional interpretation of 
Hilbert's gravitational action S f in such a way that: 



a = -GS o5 (l M /r) , l M := (Af/E ) • (10) 

and the acceleration in the full MONDian regime and the 
corresponding Tully-Fisher relation are 

1/2 

a= _ G (S M) v = G i/a s i/4 M i/4 (n) 

r 

Also, a more manageable extended fundamental quan- 
tity, directly measurable through the Tully-Fisher rela- 
tion, can be defined: 

e M := a G = 8.004 x I0" 21 m 4 s~ 4 kg" 1 , (12) 

with dimensions of velocity to the fourth over mass, for 
which 

a = ~ e -§9 (ImM , l M := (G 2 M/e ) . (13) 

With this, the acceleration of a test particle in the full 
MOND regime and the Tully-Fisher relation are: 

a = -&^, v=(e M Mf\ (14) 
r 

The choice of a new fundamental constant of nature 
has many ways in which it can be introduced into the 
theory :34]. In this work, the use of ao is kept as it is 
traditionally done, but we note the fact that eo is the best 



Sl = -\^GiJlJ f(x)V-gd\ (15) 

which slightly differs from its traditional form (see e.g. 
S HI] ) since the following dimensionless quantity has 
been introduced: 



X-=L 2 M R, (16) 

where R is Ricci's scalar and Lm defines a length fixed 
by the parameters of the theory. The explicit form of the 
length L has to be obtained once a certain known limit 
of the theory is taken, usually a non-relativistic limit. 
Note that the definition of \ gives a correct dimensional 
character to the action (flB"]) . something that is not com- 
pletely clear in all previous works dealing with a metric 
description of the gravitational field. For f(x) — X t ne 
standard Einstein-Hilbert action of general relativity is 
obtained. 

On the other hand, the matter action has its usual 
form, 



Sm = J £m\/ 3 gd 4 :r, (17) 

with £ m the matter Lagrangian density of the sys- 
tem. The null variations of the complete action, i.e. 
5 (Sh + S m ) = 0, yield the following field equations: 
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fix) X»» ~ g f{x)g ^ ~ L " M ~ 9 ^ A) f ' {x) 



irrGL 2 



4 J-iiVi 



(18) 



where the dimensionless Ricci tensor Xuv is given by: 



(19) 



and R^v is the standard Ricci tensor. 



The Laplace- 

Bcltrami operator has been written as A := V Q V a 
and the prime denotes derivative with respect to its 
argument. The energy- momentum tensor is de- 
fined through the following standard relation: SS m = 
— (l/2c) T a p 5g a P . In here and in what follows, we choose 
a (+, — , — , — ) signature for the metric and use Ein- 
stein's summation convention over repeated indices. 
The trace of equation (fT8|) is: 



/'(x)x-2/(x) + 3L^A/'( X ) = 



lirGL 2 M 



T, (20) 



R c its radius of curvature and so, relation (|23j) essentially 
means that 



R c > r. 



(24) 



In other words, the second term on the left-hand side 
of equation (|2"2")l dominates the first one when the local 
radius of curvature of space is much grater than the char- 
acteristic length r. This should occur in the weak-field 
regime, where MONDian effects are expected. For a met- 
ric description of gravity, this limit must correspond to 
the relativistic regime of MOND. 

Under assumption (|24| , equation (|22|) takes the follow- 
ing form: 



iwGM 



3bc 2 rL 



2(6-1) ' 
M 



(25) 



We now recall the well known relation followed by the 
Ricci scalar at second order of approximation at the non- 
relativistic level (19| : 



where T :— T™. 

In order to search for a MONDian solution, Bernal 
et al. 5| analysed the problem in two ways. First by 
performing an order of magnitude approach to the prob- 
lem, and second, by doing a full perturbation analysis. 
Since the second technique is merely to fix constants of 
proportionality of the problem, their order of magnitude 
approach and its consequences are discussed in the re- 
main of this section. Also, since we are interested at 
the moment on a point mass distribution generating a 
stationary spherically symmetric space-time, the trace 
equation (|20[) contains all the relevant information relat- 
ing the field equations. At this point it is also useful to 
assume a power law form for the function 



fix) = x b 



(21) 



An order of magnitude approach to the problem means 
that d/dx ~ A « — 1/r 2 and the mass density p « 
AI/r 3 . With this, the trace (|2"0")) takes the following form: 



X b (& - 2) - 3bL 



2 X 



(b-D 8 



ttGML 2 m 



M ' 



(22) 



Note that the second term on the left-hand side of 
equation (|22[) is much greater than the first term when 
the following condition is satisfied: 



Rr 2 < 



3b 
2-6' 



(23) 



At the same order of approximation, Ricci's scalar R w 
k = R~ 2 , where k is the Gaussian curvature of space and 



R = 



2 „ 

-V-a, 

c z 



(26) 



where the negative gradients of the gravitational poten- 
tial 4> provide the acceleration a := — felt by a test 
particle on a non-relativistic gravitational field. At order 
of magnitude, equation (|26| can be approximated as 



R 



2(b 



2a 



(27) 



Substitution of this last equation on relation (1251) gives 



t-r I SnGU V h 
3bc 2 r J 

(26-4)/(6-l) r (6-2)/(6- 



9T 2 



] Lm 2 (GM) 



1/(6-1) 



(28) 



This last equation converges to a MOND-like accelera- 
tion a oc 1/r if b — 2 = — (b — 1), i.e. when b = 3/2. Also, 
at the lowest order of approximation, in the extreme non- 
relativistic limit, the velocity of light c should not appear 
on equation (f2"5|) and so, the only way this condition is 
fulfilled is that Lm depends on a power of c, i.e. 



L- 2 oc 



-2b)/(6-l) 



, and so, Lm oc c . 

(29) 

As discussed by Bernal et al. @, the length Lm must 
be constructed by fundamental parameters describing the 
theory of gravity and since the only two characteristic 
lengths of the problem are the mass-length Im and the 
gravitational radius 
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GM 



(30) 



then the correct dimensional form of the length Lm is 
given by 



L M = Qr^l M , with a + p = l, 



(31) 



where the constant of proportionality £ is a dimensionless 
number that can be found by a full perturbation analysis 
technique and is given by [5[ : 



c 



2V2 



(32) 



Substituting equation pip and the value 6 = 3/2 into 
relation (|2l?]) . it then follows that 



a = = 1/2, 



i.e. 



„l/2,l/2 



(33) 



If we now substitute this last result and the value b = 
3/2 in equation (f28| we get: 



1.5 - 



0.5 



0.5 



1.5 



FIG. 2. The figure shows the transition function f(x), as 
a function of the dimensionless Ricci scalar \, for different 
regimes of gravity, converging to f(x) = X f° r X ^> 1 (general 
relativity) and to f(x) = X 3//2 f° r X ^ 1 ( a relativistic regime 
with MOND as its weak field limit) -see equation (|37p . The 
thick dash-dot curve is the extreme limiting value p — > 00, i.e. 

fix) = X 3 ^ 2 f° r X < 1 an( i /(x) — X f° r X > 1- The curves 
above and below this extreme function represent values of 
p = 3, 2, 1, for the minus and plus signs of equation (FT} 
respectively. The extreme limiting curve has a kink at x = 1- 



(a GM) 1/2 



(34) 



which is the traditional form of MOND for a point mass 
source (see e.g. [1, H?], and references therein). Also, 
the results of equation (jM)) in (|2T|) mean that 



R 



and so, inequality (|24[) is equivalent to 



Zm > r K 



(35) 



(36) 



The regime imposed by equation (j3"6")l is precisely the one 
for which MONDian effects should appear in a relativis- 
tic theory of gravity. This is an expected generalisation 
of the results presented in section [TTJ Note that in the 
weak held limit regime for which Im <C t together with 
equation (j3"6"|) yields r ^> Im ^S> r g . In this connection, 
we also note that Newton's theory of gravity is recovered 
in the limit Im 3> r r g . 

In exactly the same way as it was done to build the 
transition function for the case of extended Newtonian 
gravity in section UH a general function f(x) can be con- 
structed: 



/(x) = x 3/2 



1 ± x p+1 

I ± x 3/2+p 




forx«l, (3?) 
for x > 1- 



In other words, general relativity is recovered when \ 3* 
1 in the strong field regime and the relativistic version of 
MOND with x 3 ^ 2 is recovered for the weak field regime 
of gravity when x 1 (see Figure [5]). The unknown pa- 
rameter p > — 1 needs to be calibrated with astronomical 
observations, in an analogous form as the calibration of 
the parameter n in equation ([7]) was done. This is a much 
harder task and a matter of future research. However, 
since the non-relativistic approach to gravity explained 
in section [TT| means that the transition from the Newto- 
nian to the MONDian regimes of gravity is very sharp, 
it most probably means that the function f(x) = X f° r 
X > 1 and that /(%) = x 3 ^ 2 f° r X ^ 1) but this has to be 
tested by some astronomical observations. 

The mass dependence of \ an d Lm mean that Hilbcrt's 
action (fTS]) is a function of the mass M . This is usu- 
ally not assumed, since that action is thought to be 
purely a function of the geometry of space-time due to 
the presence of mass and energy sources. However, it 
was Sobouti [35| who first encountered this peculiarity 
in the Hilbert action when dealing with a metric gen- 
eralisation of MOND. Following the remarks by Sobouti 
[35| and Mendoza and Rosas-Guevara jH} one should not 
be surprised if some of the commonly accepted notions, 
even at the fundamental level of the action, require gen- 
eralisations and re-thinking. An extended metric theory 
of gravity goes beyond the traditional general relativity 
ideas and in this way, we need to change our standard 
view of its fundamental principles. 
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IV. F(R,T) CONNECTION 

For the description of gravity shown in section IIHI it 
follows that an adequate way of writing up the gravita- 
tional field's action is given by: 



R 



d_ 

OR' 



and 



d_ 
df' 



(43) 



The tensor 0^,, is such that Q^Sg' 11 ' :— g a ^8T a g and for 
the case of an ideal fluid it can be written as 1141: 



St = 



16ttGJ L 2 M 



(38) 



(44) 



The function Lm is a function of the mass of the sys- 
tem and in general terms it is a function of the space- 
time coordinates. For the particular case of a spherically 
symmetric space-time it coincides with the mass of the 
central object generating the gravitational field as ex- 
pressed in equations (|31~T) and (|33|) . Generally speaking 
what the meaning of M would be for a particular dis- 
tribution of mass and energy needs further research, be- 
yond the scope of this work. Nevertheless one expects 
that for a systems with high degree of symmetry (such 
as as a spherically symmetric spacetime or a Friedmann- 
Lemaitre- Robertson- Walker one) , the function M would 
be given by the standard mass-energy relation [see e.g. 
HI: 



Note that equation (|4"Tj) or (|4"2"| converge to the 
field (TT8|) and trace ((20)) relations as discussed in sec- 
tion Mil when one considers a point mass generating 
the gravitational field, i.e. when Lm = const, and so 

d/dR = L 2 M d/d x . 

In general terms, the F(R, T) theory described by 
Harko et al. [H| produces non-geodesic motion of test 
particles since: 



/ \ (SirG ( f 



M / T 



T 2 

M / T 



(V + e M „) V^ 1 In ( ) (i?,T) + V^9 Mi , 



J M / T 



(45) 



M 



-gd 3 x. 



(39) and as such the geodesic equation has a force term: 



The field equations produced by the null variations of 
the addition of the field's action Sf + S m can be con- 
structed in the following form. Harko et al. [lij have 
built an F(R 7 T) theory of gravity, so making the natu- 
ral identification: 



d 2 x^ 
ds 2 

where the four-force 



(46) 



F(R,T) := 



fix) 

T 2 

M 



(40) 



8nG , o 



1 


'8ttG 






I <? + 


(£)] 

\ M J T- 



(r-nv)v„ Pl 



(47) 



it is possible to use all their results for our particular case 
expressed in equation (|40|) . For example, the null vari- 
ations of the complete action Sf + S m for the particular 
case of equation (|40j) is given by Harko et al. [14| : 



j 2 I JX t w o T 2 J "l*" + 



J M 



2L U 
8tt(7 



T — 



and its trace is given by: 



g^uA - A^Aj, 



f_R 

T 2 

M 



(41) 



is perpendicular to the four velocity dx a /ds. As ex- 
plained by Harko et al. |14(, the motion of test particles 
is geodesic, i.e. A M = and/or V Q T Q( 3 = 0, (i) for the 
case of a pressureless p = (dust) fluid and (ii) for the 
cases in which Ft(R, T) = 0. 

In what follows we will see how all the previous ideas 
can be applied to a Friedmann-Lemaitre-Robertson- 
Walker dust universe and so, the divergence of the energy 
momentum tensor in equation (|45p is null. It is worth 
noting that this condition on the energy-momentum ten- 
sor for many applications needs to be zero, including ap- 
plications to the universe at any epoch. 



f R R V 



J 2 

M 



T 2 

M 



8ttG, 



3A(-^M=^T- [js-] T + Q 

V M J T - 

(42) 

where the subscripts R and T stand for the partial deriva- 
tives with respect to those quantities, i.e. 



V. COSMOLOGICAL APPLICATIONS 

There are many good and interesting attempts to ex- 
plain many cosmological observations using modified the- 
ories of gravity [see e.g. l3lT and references therein], how- 
ever these theories are not generally or fully consistent 
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with the gravitational anomalies shown at galactic and 
extragalactic scales discussed in sections [TT| and IIIII To 
see whether the gravitational f(\) theory developed in 
the previous sections can deal with cosmological data, 
let us now apply the results obtained in those sections 
to an isotropic Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe following the procedures first explored 
by Carranza, Mendoza, and Torres In this case, the 
interval ds is given by [2l| : 



then 



8irGF R 



l (KF F) 3HdFR 
-(RF R -F)-— — 



(53) 



where H := a(t)/a(t) represents Hubble's constant. 

With the above definitions and using the 00 component 
of the field's equations (|4"9"]l and the relation [cf . KLllj : 



ds J 



c 2 dt 2 -a 2 {t) 



dr 2 



1 — nr 2 



r 2 dn 2 



(48) 



where a(t) is the scale factor of the universe normalised 
to unity, i.e. do = 1, at the present epoch to, and the 
angular displacement dfl 2 :— d6 2 + sin 2 9 dip 2 for the po- 
lar d6 and azimuthal d(p angular displacements with a 
comoving coordinate distance r. In what follows we as- 
sume a null space curvature k = at the present epoch in 
accordance with observations and deal with the expan- 
sion of the universe dictated by the field equations (|4T|). 
avoiding any form of dark unknown component. Since we 
are interested on the compatibility of this cosmological 
model with SNIa observations, in what follows we assume 
a dust p = model for which the covariant divergence 
of the energy- momentum tensor vanishes, and so as dis- 
cussed in section HVl the trajectories of test particles are 
geodesic. 

To begin with, let us rewrite the field equations (|4Tj) 
inspired by the approach first introduced by Capozziello 
and Fang [8j (see also Capozziello and Faraoni [9]) as 
follows: 



G ^ v 



8irG 



1 + ^g Ft )Y r - + t ^ 



(49) 



where the Einstein tensor is given by its usual form: 



(50) 



and 



Ticurv . _ 

'~ 8irGF K 



(F - RF R ) - AF R 9llv 



(51) 



represents the "energy-momentum" curvature tensor. 
Since Too = pc 2 , then it will be useful the identifica- 
tion Too : = PcurvC 2 . With this last definition and using 
the fact that the Laplace-Beltrami operator applied to a 
scalar field tfj is given by [see e.g. fl9l] : 



«=4 



KC 

~7a 



(54) 



between Ricci's scalar and the derivatives of the scale fac- 
tor for a FLRW universe, then the dynamical Friedman's- 
like equation for a dust flat universe is: 



H 2 



8nG 



1 



— — Ft] -f— 

8ttG J F R 



Pc 



(55) 



The energy conservation equation is given by the null 
divergence of the energy-momentum tensor: 



dF T 
'~dT 



(56) 



For completeness, we write down the correspondent 
generalisation of Raychadhuri's equation for a dust flat 
universe: 



2- + H 2 

a 



8irGp c 



where the "curvature-pressure" 



Pcurv • — ^ c Pcurvj 



and 



w 



(F - RF R ) (2 + d 2 F R /dt 2 + 3HdF R /dt 
c 2 (RF R -F)/2- 3HdF R /dt 



(57) 



(58) 



(59) 



On the other hand, note that the mass M that appears 
on the length Lm must be the causally connected mass at 
a certain cosmic time t, since particles beyond Hubble's 
(or particle) horizon with respect to a given fundamental 
observer do not have any gravitational influence on him. 
In other words: 



M 



-9 d 3 x 



4 c 
—np 



3 



(60) 



where 



Atp = 



1 



(52) 



H(ty 



(61) 



8 



is the Hubble radius or the distance of causal contact at 
a particular cosmic epoch [2lJ . From this last relation it 
follows that the length (|3"Tj) is given by: 



(0) 



3H 2 p 
8ttG ' 



(71) 



J->M — 



ca 



1/4 #9/4 : 



In other words, the value of Milgrom's acceleration 
(62) constant a at the current cosmic epoch is such that 



and so, by using relation ([2lj) and the standard power-law 
assumptions: 



a(t) = a(t ) - 
to 



pit) = Po 



a(t ) 



(63) 



for the unknown constant powers a and /3, it follows that: 



dF K 
dt 



= b(b-l)R°- L L 



b -l r 2(6-1) 



A I 



H 



j-q-2 3/ 3 
1 — q 2 V a 



6 r 2b-2 



dFr = 3 _ s R L M 
dt 2 l ' pc 2 



(64) 
(65) 



where 



q(t) r- 



adt 2 



R- 



and 



. 1 d 3 a 
a at-* 



(66) 

are the deceleration parameter and the jerk respectively. 

With these and the value of Lm from equation (1S"2|) . 
the curvature density (|53j) is given by: 



3H 2 
8irG 



(b-l) 



(1-9) 



j - q 



2 3/3 



Substitution of the previous relations on Friedmann's 
equation ([55)1 gives: 



a « c x i? . 



(72) 



The numerical coincidence between the value of Mil- 
grom's acceleration constant ao and the multiplication 
of the speed of light c by the current value of Hubble's 
constant Hq has been noted since the early development 
of MOND [see e.g.[l3|, and references therein]. Note that 
equation (|72p means that this coincidence relation occurs 
at approximately the present cosmic epoch in complete 
agreement with the results by Bernal et al. 0] where it is 
shown that ao shows no cosmological evolution and hence 
it can be postulated as a fundamental constant of nature. 

For the power law (|21[) and the assumptions made 
above, it follows that the energy conservation equa- 
tion (f5"6")) is given by: 



(p + 3Hp) 



where: 



8ttG 



A- 



RH 



(73) 



R 



9 b- 
2~b 



A :-- 
27 (6 



9 

if 



3 6(6-l)(j-g-2) 
2 



Direct substitution of the density power law (l€>3[) into 
relation (|73|) gives a constraint equation between a, /3 
and 6: 



ff2 = SvrGp 



(68) 



where 



l + (b-l) 



j-q-2 4(l-g) 3/3 
1 - g 6 2 ^ P a 



(69) 

is a dimensionless function. 

An important result can be obtained evaluating equa- 
tion at the present epoch, yielding: 



a = 



( 36 _5)/(b_l) 



ci? 
J (70) 

where the density parameter f^itt & t the present epoch 
has been defined by it's usual relation: 



1 /9- 
a V 36 



56 



(74) 



Let us now proceed to fix the so far unknown param- 
eters of the theory a, (3 and 6. To do so, we need reli- 
able observational data and as such, we use the redshift- 
magnitude SNIa data obtained by Riess et al. [33[ and the 
following well known standard cosmological relations [see 
e.g. Hl|: 



1 + z = a(t )/a(t), (75) 
li (z) = 5 log 10 [H a d L (z)} - 5 log 10 h + 42.38, (76) 

c 



d L (z) = (1 + z) 



o H(z) 



dz, 



(77) 



for the cosmological redshift z, the distance modulus p, 
the luminosity distance and where the normalised 
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FIG. 3. Redshift magnitude plot for SNIa showing the dis- 
tance modulus fi as a function of the redshift z for SNIa as 
presented by Riess et al. [3^1 - The dotted red line shows the 
best fit to the data with the /(x) gravity theory applied to a 
flat dust FLRW universe (see text) with no dark components. 
The continuous blue line represents the best fit according to 
the standard concordance dust ACDM model. 



Hubble constant h at the present epoch is given by 
h := H / (lOOkms-yMpc). Also, from equation §5$) it 
follows that 



H(a) = H 



a(to) 



-l/a 



H (l + z) 1/a , 



(78) 



and the substitution of this into equation (f77|) gives the 
distance modulus as a function of the redshift z. This 
means that the redshift magnitude relation (|To| is a func- 
tion that depends on the values of the current Hubble 
constant Hq and the value of a. Figure [3] shows the 
best fit to the redshift magnitude relation of SNIa ob- 
served by Riess et al. |33[, yielding a = 1.359 ±0.139 and 
h = 0.64±0.009. The best fit presented on the figure was 
obtained using the Marquardt-Levenberg fit provided by 
gnuplot (http://www.gnuplot.info) for non-linear func- 
tions. These values do not provide the whole description 
of the problem, since /3 and b are still unknown. How- 
ever, according to the constraint equation (|74|) only one 
of them is needed in order to know the other once a is 
known. 

The parameter f3 can be found from conservation of 
mass arguments, since the total mass of the universe 
M iot — An J flmax pr 2 a(t) 3 dr = const., where the upper 
limit of the integral is the radius of the whole universe. 
Since a(t) and p{t) are time dependent functions, the 
only way the mass of the universe is conserved is by re- 
quiring a 3 p — const, and so, ft = 3. This argument is 
exactly the one used in standard cosmology when deal- 
ing with a dust FLRW universe [see e.g. [21]. Using this 
value of /3 and the one already found for a, it follows 
that b = 1.57 ± 0.56, which is within the expected value 



of b = 3/2 discussed in section HIT1 

For completeness, we write down a few of the cosmog- 
raphycal parameters obtained by this /(x) gravity ap- 
plied to the universe: 



(79) 



h = 0.64 ±0.009, q = -0.2642 ± 0.075, 
jo = -0.1246 ±0.004. 

VI. DISCUSSION 



As explained by Carranza, Mendoza, and Torres [To| . 
the obtained value b « 3/2 is a completely expected re- 
sult due to the following arguments. As explained in 
section HH a gravitational system for which its character- 
istic size r is such that x := I m / r < 1 \s in the MONDian 
gravity regime. For the case of the universe, x ~ a few 
and as such if not totally in the MONDian regime of 
gravity, then it is far away from the regime of Newtonian 
gravity. The relativistic version of this means that the 
universe is close to the regime for which /(x) = X 3 ^ 2 and 
so b = 3/2. This is a very important result since seen 
in this way, the accelerated expansion of the universe is 
due to an extended gravity theory deviating from gen- 
eral relativity. It is quite interesting to note that the 
function fix) = X 3 ^ 2 which at its non-relativistic limit 
is capable of predicting the correct dynamical behaviour 
of many astrophysical phenomena, is also able to explain 
the behaviour of the current accelerated expansion of the 
universe. 

Seen in this way, the behaviour of gravity towards the 
past (for sufficiently large redshifts z) will differ from 
fix) = X 3/2 and eventually converge to fix) = X, i-e- 
the gravitational regime of gravity is general relativity for 
sufficiently large redshifts. A very detailed investigation 
into this needs to be done at different levels in order to 
be coherent many different cosmological observations [see 
e.g.[20[. This in turn can serve to calibrate the index p of 
the transfer function /(x) as presented in equation ([57]) . 
which has a very soft transition when p = —1, i.e., 



fix) 



^3/2 



1 + x 3/2 ' 



(80) 



and also has a very sharp transition when p —¥ oo, with 
the step function: 



fix) 



( 3 / 2 , for < x < 1, 
6 for X > 1- 



(81) 



In this respect, perhaps something close to a sharp tran- 
sition (|81l) will be observed since, as mentioned in sec- 
tion|TTl at the non-relativistic level different astrophysical 
observations show a sharp transition from the Newtonian 
to the MONDian regimes. This sort of decision has to 
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be taken with care and such a full description requires to 
analyse in full detail the whole Friedmann-like equations: 



hrG 



F 7 



dF T 



1 / dF T 



U 2 SvrG 




H = 3 




a n kc 2 


2-+H^ + - 




a 


a 2 



dt 

c 4 F T \ p 
8nG ) F R + Pcur 
8ttG P 2pc 2 F T 



TP dp 



(82) 



^r, (83) 



8nG Pc 



Fb 



(84) 



These equations are directly obtained from taking the 
null covariant divergence of the energy momentum ten- 
sor, the 00 component of the field equations (|49l) and the 
density p contains all species of matter and/or radiation. 
The curvature density /9 C urv and the curvature pressure 



Pcurv are related to one another by relation (|58[) with ui 
given by equation (|5§)) . 

It is quite remarkable that a metric extended theory 
of gravity is able to reproduce phenomena from mass 
and length scales associated to the solar system up to 
cosmological scales. There are many more astrophysical 
challenges that this theory needs to address, in particular 
with respect to lensing at different scales and the dynam- 
ics associated to galaxy clusters. These will be addressed 
elsewhere. 
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